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Brown dwarfs and giant gas planets are substellar objects whose spectral appearance is de- 
ON \ termined by the chemical composition of the gas and the solids/liquids in the atmosphere. 

Atmospheres of substellar objects possess two major scale regimes: large-scale convective 
motions + gravitational settling and small-scale turbulence + dust formation. Turbulence 
initiates dust formation spot-like on small scale, while the dust feeds back into the turbulent 
fluid field by its strong radiative cooling. Small, imploding dust containing areas result which 
eventually become isothermal. Multi-dimensional simulations show that these small-scale 
dust structures gather into large-scale structures, suggesting the formation of clouds made of 
dirty dust grains. The chemical composition of the grains, and thereby the chemical evolution 
of the gas phase, is a function of temperature and depends on the grain's history. 



in 
o 



Of 

6 

CO 

C3 



X 

5-S 



1 Introduction 



The first brown dwarf Gliese 229B has been discovered 10 years ago by direct imaging (Kulka- 
rni & Golimovsky 1995). These faint (L* = 1(T 7 . . . 10 _1 L Q ), cool (T cfT = 500 ... 3200 K), 
and small (M* = 0.01 . . . O.O8M ) objects bridge the physical and chemical gap between the 
classical understanding of stars (M* > O.O8M ) and planets (M* < O.OIMq). Much closer 
by, direct images revels spotty, cloudy, and vortex surface pattern in the giant planet atmo- 
spheres in our own solar system (by the Cassini and Galileo spacecrafts fly-bys of Saturn and 
Jupiter, respectively) which guide our imagination for substellar but extra-solar atmospheres. 
The other major source of information is the measurement of the energy distribution of the 
stellar radiative flux emerging from the object's atmosphere (e.g. for Gliese 229B Oppen- 
heimer et al. 1998). The interpretation of the resulting spectral energy distribution demands 
a certain compleatness 1 of the adopted substellar atmosphere model. Substellar atmospheres, 
i.e. giant gas planets and brown dwarfs, are very cool and therefore exhibit a rich molecular- 
and solid-/liquid-phase chemistry. Transitions between the phase regimes are to be expected. 
Therefore, models of substellar atmospheres - as the interface to the physical and the chemical 
state of the object - need to represent the circuit of dust (Helling 1999, Woitke & Helling 2003) 
which includes the formation of dust, the chemical composition of dust and gas, gravitational 
settling (rain), its feedback on the dust formation process, and element replenishment by 
upward convective motions in addition to hydrodynamics and radiative transfer. In contrast 

1 e.g. in modeling the molecular regime, chemistry and hydrodynamics 
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Table 1: Time scales of the processes involved in the circuit of dust in a substellar atmosphere. 
Temperature T, density p, pressure scale height H p , convective velocity f C onv> an d velocity of 
sound c s are adopted from the model results by Allard et al. (2001) and Tsuji (2002). For 
more details on r s i n k see (Woitke & Helling 2003). 



to terrestrial planets which possess their solid surface as continues source of seed particles 2 , 
the actual formation of the first (solid or liquid) surface out of the gas phase has to be con- 
sidered in substellar atmospheres. Convection is an efficient mechanism to continuously and 
intermittently dredge up fresh, uncondensed gaseous material from the very bottom of the 
atmosphere. The convection furthermore serves as turbulence engine inside the atmosphere. 
Consequently, modeling and understanding a substellar atmosphere means to model and to 
understand a reactive, dust forming, turbulent fluid field. 

2 Catching the small scales 

Classical models for substellar atmospheres represent the whole turbulent scale spectrum by 
only one scale, the mixing length. These models have given very reasonable fits to observed 
spectra in certain wavelength regimes but are challenged by the progress in observational 
techniques which lead to observations e.g. with higher resolutions and at longer wavelength 
(A > 12/um). Other models like Reynolds stress and LES are in progress, all of them being 
challenged by the closure problem, i.e. the treatment of the smallest, unresolved scales. 

In order to provide insight and understanding of the small scale regimes of a substellar 
atmosphere, the interaction of turbulence and dust formation has been studied by utilizing 
ID and 2D simulations in the present work. The general phenomenology of a substellar at- 
mosphere model can be demonstrated by estimating the time scale of the individual processes 
(Table 1). The gravitational settling time scale of grains r s i n k is the largest and is compara- 

2 Seed particles on Earth are called aerosols which are for instance volcanic dust and tire particles, or they 
come from fire in the tropics and from smoke-tracks. 
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ble to a typical convective mixing time for large grains. The smallest time is needed by the 
formation of seed particles out of the gas phase. The dust growth time scale is of the order 
of the crossing time of an acoustic wave. Hence, two time regimes appear: (i) a quasi-static 
regime governed by gravitational settling and large-scale convective motions, and (ii) a dy- 
namic regime governed by the dust formation and small-scale waves (turbulence). Actually, 
a third regime (hi) is to be faced which concerns the computing time needed to resolve the 
turbulent problem (last entry Table 1). In order to tackle regime (iii), the small-scale regime 
(ii) was investigated. Here, gravitational settling can be neglected and convection acts only 
indirectly as turbulence driver. 

The following system of dimensionless equations has been solved where Eqs. 1- 3 are the 
equation of continuity, of motion, and the energy equation, respectively. The source term in 
Eqs. 3 is due to radiative cooling modeled by a relaxation ansatz. 



(p)t + V-(pv) = (1) 

1 M 2 

(H( + V-(pvov) = -—VP- 7 — pg (2) 

(pe) t + V.(v[pe + P]) = Rd k (7r E — T 4 ) (3) 



; v nct 

(pLj) t + V-(vpLj) = DarSejJ, + Daf rJ 4— pLj_! (4) 

R 

(pe x )t + V • (v pe x ) = -J2« UC E1 Da r v^iVi J* + vf El Daf n x , rVrcl , x a r pL 2 ) (5) 

r=l 

Equations 4 (j = 0, 1, 2, 3; pLj(x, t) = f-^? /(V, x, t) V^ 3 dV dust moments, f(V, x, t) grain size 
distribution function) model the dust formation as two step-process, namely, seed formation 
and mantle growth/evaporation being the first and the second source term, respectively. 
Equations 5 are element conservation equation for each chemical element e x (x= Mg, Si, O, 
Fe, Al, Ti) involved in the dust formation processes, hence each source term in Eqs. 4 is a 
sink for Eqs. 5 (Helling et al. 2001 for details). A strong coupling exists between Eqs. 1- 3 (5 
equations) and Eqs. 4, 5 (11 equations) due to the dust opacity k since it changes by order 
of magnitudes if dust forms. 



2.1 Turbulence <-> Dust formation 

The feedback between turbulence and dust formation can be studied in detail by ID sim- 
ulations. Interacting turbulence elements are modeled as superimposing expansion waves. 
Figure 2 demonstrates the time evolution of the system at the site of constructive wave in- 
teraction. 

Turbulence — > Dust formation: 

At about the time of superposition (black cross on T-curve, l.h.s. 2nd row), the nucleation 
threshold temperature (here for Ti(>2 seed formation) is crossed and dust nucleation is ini- 
tiated, hence the nucleation rate J* increases. Many solid compounds are already thermally 
stable at such low temperatures which results in a very rapid growth of a mantle on the sur- 
face of the seed particles (here Mg2Si04[ s ], Si02[ s ], Fe[ s ], Al203[ s ], Ti02[ s ])- As the amount of 
dust formed increases, the opacity k of the dust-gas mixture increases by order of magnitudes. 
Therefore, the radiative cooling causes the temperature T to drop considerably (l.h.s., 2nd 
row). A classical instability establishes where the reason supports the cause. The cooler the 
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Figure 1: The feedback between turbulence and dust formation in a ID simulation of inter- 
acting turbulence elements. The time evolution of all quantities is depicted at the site of 
maximum superposition of the turbulence elements. 

Parameter: T ref = 1900K, T RE = 1634K, p ref = 10" 4 g/cm 3 , M=0.1 (tt rof = 3 • 10 4 cm/s), 
i re f = 3s, Z re f = 10 5 cm. [numerical parameter: N x = 500, 2V& = 500, Ax = 3.94 • 10 2 cm, 

^max — Z re f/2] 

1st row: l.h.s. - p density, r.h.s. - p pressure; 2nd row: l.h.s. - T temperature (cross 

- time of maximum wave superposition), r.h.s. - u fluid velocity; 3rd row: l.h.s. 
log J*/n<H> nucleation rate [1/s], r.h.s. - log (a) mean grain size [cm]; 4th row: l.h.s. - 
Vtot = E V B cumulative volumes {%] (orange - % g si0 3[8] , brown - 14i g Si0 3[s] + %0 2H , green 

- Hi g Si0 3[s] + ^Si0 2[s] + ^Fe [s] , light blue - ^M gS i0 3[a] + V S io 2[e] + V Fe[s] + Vai 2 o 3[s] , dark blue - 
^M g Si0 3[s] + Vsi0 2[s] + ^Fe [s] + Vai 2 3[b] + %i0 2[s] ), r.h.s. - log n y number density of gaseous key 
species for dust formation [1/cm 3 ] (orange - Mg, brown - SiO, green - Fe, light blue - AlOH, 
black - TiO a ). 
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gas, the more dust forms, the faster the temperature drops etc. This run-away effect stops 
if all condensible material was consumed or if the temperature is too low for further efficient 
nucleation. The time of maximum nucleation rate corresponds to a minimum in mean grain 
size because suddenly the available gaseous material is needed for a much larger number of 
grains. 

Dust formation — > Turbulence: 

The strong temperature gradient causes a strong raise in density by a moderate pressure 
gradient. Without such a strong cooling, pressure equilibrium should adjust. Hence, the dust 
forming areas implode and cause a considerable disturbance of the velocity field (here up to 
10%). A feedback-loop turbulence =>■ dust formation =>■ turbulence established as result of 
non-linear coupling of the model equations. 
Chemistry: 

Figure 2 depicts in the lowest two panels the strong feedback of the turbulent dust formation 
process on the chemical composition of the dust grains and on the remaining gas phase. The 
cumulative dust volumes (l.h.s., 4th row) show that the dust composition changes according to 
the temperature from Mg 2 Si04[ s ]/Si02[ s ]/Fe [s] -rich (40%/40%/10%) with Al and Ti - oxide 
impurities to a mean grain composition of 55% SiC"2[ s ], 25% Fe[ s ], and 10% Mg2Si04[ s ] with 
Al and Ti - oxide impurities. The chemical composition of the dust is also imprinted in the 
element abundances of the gas phase and the molecular abundances adjust accordingly. 

2.2 Dust clouds growing from small — > lage 

ID simulations allow a very detailed investigation of the physio-chemical interactions but 
fail to explain 2- or 3-dimensional phenomena like cloud or vortex formation as it was sug- 
gested by observations of Saturn's and Jupiter's atmospheres. The study of the formation and 
the possible appearance of large scale dust clouds in substellar atmospheres was performed 
utilizing a model for driven turbulence. Turbulence is modeled by the superposition of Nk 
modes each having a Kolmogoroff velocity amplitude (for details Helling et al. 2004). Con- 
vection is believed to drive the turbulence in a real substellar atmosphere. The 2D simulation 
(Fig. 2.2) is started from a homogeneous, dust free medium which constantly is disturbed by 
the turbulence driving from the left, the right and the bottom side during the simulation. 

During the initial phase of the simulations small scale nucleation events occur where 
turbulence causes the local temperature to drop below the nucleation threshold (compare 
Sect. 2.1). Observe that a locally maximum J* in Fig. 2.2 a) is immediately followed by 
an subsequent increase in n& (panel b). These dusty areas tend to increase as the fluid 
motion transports the dust into areas with a still undepleted gas phase. And indeed, as the 
simulation proceeds in time, larger and more compact cloud-like dust structures are formed 
(panel c). These large-scale structures are the result of the hydrodynamic fluid motion which 
gathers more and more dust also by the vortices appearing in the velocity field. Note the 
mushroom- like structure evolving e.g. in the right lower corner in panel c). Strong radiative 
cooling causes these dusty areas to become isothermal. Eventually, the cloud will leave the 
test volume or it will get disrupted and the small-scale fragments move out of sight. The 
whole dust cloud formation cycle can start again only if metal species (Mg, Si, Fe, Al, Ti, O) 
are replenished from outside like it has to be expected to occur by convection in a substellar 
atmosphere (panel d). 
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Figure 2: T rcf = 2100K, T RE = 1980K, g = 10 5 g/cm 2 , /9 ref = 3.16 
(u ve { = 3 • 10 5 cm/s), t Te f = 0.3, Z re f = 10 5 cm [numerical parameter: 
N k = 500, Ax = 3.94 • 10 2 cm, Z max = Z ref /2] 

log rid number of dust particles - false color (red = max, blue = min) , log J* nucleation rate 
- contour lines (green = max, red = min), v = v x + v y fluid velocity - vector arrows. 
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3 Conclusion 

Substellar atmospheres, i.e. atmospheres of giant gas planets and brown dwarfs, possess 
two scale regimes: (i) the large-scale convection causing element replenishment of the upper 
atmosphere and the counteracting gravitational settling of dust causing an element depletion 
of the upper atmosphere 

(ii) small-scale turbulence and dust formation establishing a feedback loop. 

Regime (i) can be considered as quasi-static with view on the dust formation process. The 
dynamic regime (ii) is determined by the turbulence initiating spot-like dust formation on 
small spatial scales. The strong radiative cooling by dust results in an implosion of the dust 
forming areas which eventually become isothermal. Velocity disturbances occur which feed 
back into the turbulent fluid field. The resulting mesoscopic flow gathers the dust in even 
larger, more compact cloud-like structures. The final chemical composition of the grain is a 
function of the grains history. 
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